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Abstract 

We study time correlation functions of operators representing heavy-light K-A like systems 
at various relative distances r. The heavy quarks, one in each hadron, are treated as static. An 
anisotropic and asymmetric lattice is used with Wilson fermions. Our goal is to extract an adiabatic 
potential and thus learn about the physics of the five-quark system viewed as an hadronic molecule. 
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INTRODUCTION 



Lattice QCD studies of hadron-hadron interactions are the gateway to nuclear physics 
through first principles . From a lattice simulation point of view the nucleon-nucleon 
interaction is probably the most challenging case. This is evident given the large spatial 
size of the deuteron, and in particular, the insight that the physics of the long-range strong 
interaction is driven mostly by the pion cloud |3[ . The latter may be taken as an indication 
that chiral symmetry and a full (unquenched) lattice action are high priorities for simulations 
aiming at quantitative results. 

Aside from the above prominent case, however, interactions in other two-hadron systems 
are worth investigating as well, because this might lead to new insights into the structural 
features of some of the experimentally known baryon resonances j4|. In particular, we here 
ask if some of those may be understood as hadronic molecules, similar to the deuteron, 
but possibly with different physics of the interaction mechanisms in which quark and gluon 
degrees of freedom play a role. Prime candidates for such systems are pairs of hadrons con- 
taining one heavy quark each because, in the spirit of the Born-Oppenheimer approximation, 
the (slow) heavy quarks naturally serve to define the centers of two hadrons while the (fast) 
light quarks and gluons provide the physics of the interaction. Studies along those lines have 
been done before in the context of meson- meson and baryon-baryon systems 0,0,0]. 

We here report on the current status of interaction studies of a heavy-light meson-baryon 
(five-quark) hadron with the quantum numbers of an s-wave K-A system. Because the 
static approximation is employed for the heavy-quark propagator, the total energy can be 
computed as a function of the relative distance r between the heavy quarks. Thus an 
adiabatic (Born-Oppenheimer) potential V a {r) is extracted. This can be used to address the 
possibility of molecule-like structures. 



SIMULATION DETAILS 



Two-hadron interpolating fields are constructed from standard local operators for the K + 
and the A particles at relative distance r and projected to total momentum zero 



O a (f; t) = V- 1 ' 2 E V, x - yK + (xt)A° a (yt) 



(1) 
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FIG. 1: Diagrammatic representation of (@J). Thick (black) lines indicate heavy-quark propagators, 
and thin (blue and red) lines depict light (u and d) quark propagators. 

Here V is the spatial lattice volume and a is a Dirac spinor index. Then, with 

a M (f;t) = Ot(f;t) 74iQM (2) 



the correlation function 



(3) 



where s is the relative distance at the source, can be expressed in terms of fermion propa- 
gators. The flavor assignment K + A° ~ su uds causes the separable term to vanish. Writing 
H(xt, yt ) and G(xt, yt ) for the heavy (s) and light (u,d) quark propagators, respectively, 
one obtains 



C = (^[if(#,jH7$)ir(Fi + st ,^ 
y 

G(yt, rit )[G(yt, f x t )G{y + fit, f x + st ) - G(yt, f x + st )G(y + fit, nt )]) 



(4) 



For clarity the rather involved color and spin index structure is not shown in (£Q). Also, 
translational invariance has been used to arrive at the above expression, and an arbitrary 
space site fi was introduced in this context. A diagrammatic representation of Q is shown 
in Fig. [TJ 

The heavy-quark propagators are employed in the static approximation. For (unim- 
proved) Wilson fermions with hopping parameter k this means that the propagator is taken 
in the limit k — > , resulting in 



H(xt, yt ) = <% ? g (2k)*-* -(1 + 74) U\x] t t) , 



(5) 



where U(x; tot) is the product of SU(3) link variables along a straight line from (xto) to (xt) 



The distance r = is rather special |(| because a color singlet operator, as realized by 
can also be achieved by a "color twisted" version of (0) where quarks in K + and A are 
combined into a color singlet. Because we do not consider color twisted operators in this 
work we restrict ourselves to non-zero relative distance. Thus, because H(yt,y + ft) oc 
and H{fi + st , fit ) oc only the last two diagrams in Fig. [T] make a contribution to the 
correlation function for non-zero relative distance. By way of (JSJ) those diagrams in Fig. ^ 
are proportional to S^S^g. Thus (jlj) becomes 

C = 6 f > f 2 _f i (#(fit,fit )if(f 2 to,r2t) (6) 
G(f l t,f l t )[-G(f 1 t,f l t )G(f 2 t,f 2 t ) + G(f l t,^ for f ^ 0, 

where f 2 is yet another arbitrary space site. We observe that sources at fixed spatial sites 
only are needed. However, an undesirable consequence of the static approximation is that 
the site sum has vanished from (@J) and thus is no longer working to improve statistics. 

The final correlator we use is extended from © to a matrix by employing several levels 
k = 1 . . . K of operator smearing. The procedure amounts to replacing in (0) all light-quark 
fields ip, with smeared fields ^ k \ r ijj^ k \ The computation of light-quark propagators thus 



requires various leve 
gauge field fuzzing 



s of smearing at the source and at the sink. We have used APE-style 



9j and Wuppertal fermion smearing with common values for the 
strength parameters a = 2.5 and the number k of iterations. No smearing, nor link variable 
fuzzing, was done for the heavy, static, quark fields in order to preserve spatial locality, i.e. 
the 5 factor in (jSJ). Thus, writing O — > = 0[ip^ k \^ k ^ . . .], the correlator © becomes 
a K x K matrix 

C M (f;Mo) = (Oi fc} (r-;t)Of(f;t )) with k,£ = l---K. (7) 

The expression for C M in terms of quark propagators still has the form given by ©, however, 
light propagator elements are replaced with smeared ones, G — > G M , with appropriate 
smearing levels at source and sink. The correlation matrix ((7j) is hermitian by construction. 

The lattice geometry is chosen as L\ x L 2 x L 3 x L 4 = 8 x 8 x 32 x 16 with bare 
lattice constants a± = a 2 = 2a^ = 2a^ in the respective directions. This choice of an 
asymmetric and anisotropic lattice provides a fine mass resolution, in t-direction, and the 
same spatial resolution for the adiabatic potential, as the static sources are placed along the 
z-direction. The positions of the latter are at x = (5, 5, n, 3) with n — 1,2, 3, 4, 8, 11, 13, 17, 



with r\ t 2 = (5,5,n) in reference to ©• In this way all possible relative distances r = \r\ 
in the range r = 1 ... 16 can be obtained. Note that periodic boundary conditions, in the 
spatial directions, allow us to only utilize half of the extent of the lattice in z-direction. 

We have used the Wilson plaquette action with Wilson fermions in a quenched simulation. 
The gauge field couplings in the /i-z/ planes and the hopping parameters in directions p are 
given by, respectively, 

^ = 13 (n n \2 aIld K A* = 1^4 T " ( 8 ) 

The simulation was done at (3 — 6.2 with four values k = 0.140, 0.136, 0.132, 0.128 of 
hopping parameters using a multiple mass solver 



ANALYSIS 

The correlation matrix (J7|) is constructed with three smearing levels, K = 3. It is then 
diagonalized separately on each timeslice t, using singular value decomposition. For asymp- 
totic t the largest eigenvalues correspond to the ground state of the two-hadron system. The 
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is used to extract masses. 



maximum entropy method (MEM), as implemented in 

In order to set the physical mass scale the nucleon (N) and vector meson (p) masses are 
fit to the model y = C1 + C2X + C3 ln(l + x), with x = (am n ) 2 , using the four available hopping 
parameters. Here o = 03 = 04 is the common lattice constant in the z and t directions. The 
extrapolated values of am^ and CLTflp^ cLS X — y 0, are used to set the reduced mass of the N-p 
system to its experimental value m Np = 424MeV. This gives o = 0.096fm (a -1 = 2055MeV) 



At the time of this writing the analysis is limited to 90 gauge configurations and to the 
largest value of the hopping parameter, k = 0.140. The corresponding results for the total 
energy W are shown in Fig. |2K- The uncertainties are the (rms) widths of the spectral peaks 
emerging from the MEM analysis. Due to periodic boundary conditions the data points 
replicate for r > 16a. For a preliminary analysis we use the simple model 

aV(x) = exp(— a!X 2 )(a 2 + a 3 x 2 + a 4 x 4 ) with x — r/a. (9) 

The fit to the data is then made using the periodic replication 

aV L {x) = a Q + aV{x) + aV{L 3 -x), L 3 = 32 , (10) 
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FIG. 2: Total energy W of the heavy- light K-A like system versus the relative distance r with 
MEM based uncertainties (a). The continuous lines represent fits with Q and (|lUj) . The adiabatic 
potential V a is shown in physical units (b). Results are for k = 0.140. 

of the model, which then has five parameters . . . a^. Figure^ shows the corresponding 
result. The plot in Fig. Efc> represents Q in terms of the physical scale, the adiabatic 
potential is V a {r) = V(r/a). Pending an anticipated increase of gauge configurations no 
error analysis for V has been done at this time. However, error bands may be inferred from 
the uncertainties on the data points in Fig. |2k- 

The attractive dip of V a (r) at around r « 0.8-1.2fm, in Fig. Efc>, may be significant 
enough to produce a molecule-like structure. As a first attempt to calculate phase shifts, 
we solve a standard non-relativistic scattering problem (Schrodinger equation) employing 
the computed adiabatic potential V a with several values tur for the reduced mass. For the 
' atter we use the (experimental) reduced masses of the systems K-A, D-A c , and B-A^ from 
"\. The resulting s-wave scattering phase shifts So(p) are shown in Fig. El According to 
Levinson's theorem (o^(0) — 5e(oo) = rnc) the number n of bound states is zero for the K-A 
system. However, there is indication of an emerging resonance (rising phase shift) somewhat 
below p = 0.1a -1 , or E = 62MeV in terms of the relative kinetic energy. It should be kept 
in mind that Fig. 01 reflects a result at m n > (k = 0.140). At present the strength of this 
feature as m v — > is an unresolved question. Also, relativistic effects are not taken into 
account. The latter are less significant for the D-A c and B-Af> systems, respectively. Those 
appear to be bound by the adiabatic potential. 
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FIG. 3: Scattering phase shifts (s-wave) for the K-A, D-A c , and B—Af, systems (black, red, and 
green, respectively) versus the, nonrelativistic, relative momentum p. The arrow indicates the 
relative kinetic energy for the K-A system at p = 0.1a -1 . 

ASSESSMENT 



Although in its present state the simulation is not conclusive, the results give a hint at 
potentially interesting physics. It appears conceivable that the known hadron mass spectrum 
may contain five-quark hadrons with a molecule-like structure. Our preliminary results 
would point to a resonant state with an excitation energy typical of a nuclear system, say 
~ 50±50MeV. As a possible candidate with the appropriate quantum numbers the N(1650) 
comes to mind |4j. Its mass lies just 40MeV above the K-A threshold. However, we should 
caution that the extraction of masses typical for nuclear physics from a lattice simulation is 
difficult, because residual hadron-hadron interactions are, at least, one order of magnitude 
less that baryon rest masses. Therefore, awaiting the analysis with a reasonable number of 
gauge configurations, the above scenario should be considered an interesting possibility. 



* This material is based upon work supported by the National Science Foundation under Grant 
No. 0300065 and upon resources provided by the Lattice Hadron Physics Collaboration LHPC 
through the SciDac program of the US Department of Energy. 
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